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INTRODUCTION 


The  specific  purpose  of  this  report  i 
through  the  use  of  tables,  to  calculate  mi 
component  at  a  given  confidence  level  for 
reliability. 

It  is  assumed  that  a  representative  s 
is  available}  the  sample  size  must  be  two 
statistical  theory  underlying  this  report 
ref.  (1). 


s  to  enable  one  , 
ssion  life  of  a 
a  given  confiden 

ample  of  failure 
or  greater.  The 
is  available  in 


ced 


s 
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2.  NOTATION 

rci)  -  gamma  function  of  x 


zf-  <? 


=  standard  normal  distribution  function  of  x 


--t"- 

U  e 


fZ-TT 


^  =  standard  normal  fractile  function  of 

or  the  inverse  of  £  m\  =  o( 


n  =  sample  sixe 
X  =  mission  life 

Rc  =  Confidenced  Reliability  (lower  confidence 
bound  on  population  reliability) 

R*  =  Reliability  point  estimate 
n- 

Maximum  Likelihood  estimates  of  normal 


parameters 

Maximum  Likelihood  estimates  of  Lognormal 


parameters  (sample  mean  and  standard 
deviation  of  natural  logs  of  data  points)* 
C  =  confidence  level 
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DESCRIPTION  OF  PROBLEM 


The  prob lem  is : 

given  n,  c  and  R^, ,  solve  the 


following  equation  for^. 


The  point  estimate  of  normal  reliability  is  given  by: 


^  = 


o 


(3) 


Therefore,  once  Eqn.  (1)  is  solved  for 
mission  life  x  may  easily  be  calculated 


4?*  ’  'j! 

for  any  yu. 


and 


If  one  is  restricted  to  a  certain  reliability  R^,  and  a 
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certain  confidence  level  C,  the  mission  life  calculation  is 
solved  once  and  for  all  by  solving  Eqn.  (1)  for  a  suffici¬ 
ently  large  sequence  of  values  of  n. 

The  point  estimate  of  lognormal  reliability  is  given 


can  be  used  to  calculate  mission  life  for  the  lognormal 
model  also. 

The  method  selected  to  solve  Eqn,  (1)  is  called  the 
midpoint  or  bisection  method.  This  method  is  probably  the 
simplest  and  most  straightforward  method  of  solving  nonlin¬ 
ear  equations  in  one  unknown.  The  method  consists  of: 

1.  Locating  the  root  r  between  and  X2  (by  obser¬ 
ving  a  change  of  sign  in  the  function). 
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Evaluating  the  function  at  the  middle  of  the  in¬ 


terval  (xi*  *2^  then  observing  on  which  side  of  the 

midpoint  the  function  changes  sign. 


3.  Redefining  either  or  X2  accordingly  and  re¬ 


turning  to  step  2  if  the  interval  size  is  not  sufficiently 
small. 

The  point  at  which  this  process  is  stopped  is  deter¬ 
mined  by  the  fact  that  an  easily  calculated  upper  bound  on 
the  relative  error  in  the  root  is  always  available. 


Figure  1 


Referring  to  Fig.  (1),  suppose  the  process  is  to  be  stopped 
with  the  present  values  of  x^^  and  X2.  The  final  estimate 
of  the  true  root  r  will  be  r*  =  (x^  +  X2^^^»  therefore: 
absolute  error  =  I  |f"  — 


r  —  r 
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If  the  process  is  continued  until  p*  is  less  than 
sired  fractional  error,  p  will  also  be  less  than 
and  the  process  may  be  halted.  This  is  all  true 
that  can  be  evaluated  accurately. 


some  de- 

this  error 
assuming 


The  next  problem  to  be  considered  is  evaluating  the 
integrand  in  Eqn  (1).  Formulas  26.2.12  and  26.2.17  of  ref. 

(2)  are  used  to  evaluate  the  normal  distribution  function 

i. 

Part  of  the  integrand  — 


as  n  becomes  large.  Fig.  (2)  indicates  this  behavior. 
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This  behavior  at  s*  behooves  one  to  see  that  the  quadra¬ 
ture  formula  used  to  perform  the  integration  evaluates 
the  integrand  at  s*.  A  three  point  Gaussian  quadrature 
formula  is  used  in  a  composite  manner  to  evaluate  the  in¬ 
tegral  in  Eqn.  (1).  The  principle  of  successive  consis¬ 
tency  is  used  to  examine  the  accuracy  of  the  solution  to 

Eqn.  (1);  i.e.,  the  integral  is  evaluated  with  a  succes¬ 
sively  smaller  interval  si*e  until  the  solutions  to  equa¬ 
tion  (1)  become  consistent.  The  only  remaining  aspect 
of  the  problem  is  the  evaluation  of  the  constant  in  equa¬ 
tion  ( 1) . 
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(7) 


o-c-yrf-^) 


if  n  is  even, 


(  8) 


(9) 


Expressions  (8)  and  (9)  are  easily  derived  using  the  re¬ 
cursion  formula  for  the  gamma  function  : 
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4.  DESCRIPTION  OF  COMPUTER  PROGRAM 
A.  Main 

Significant  variable  names: 

N  =  sample  sire 

NR  =  no.  of  confidenced  reliabilities 
NIT  =  no.  of  iterations  for  obtaining  consistency 
in  results 

R  =  RC  =  confidenced  reliability 

NINT  =  no.  of  intervals  between  xero  and  s*  for  com¬ 
posite  integration 

C  =  confidence  level 

FRE  =  fractional  error  tolerated  in  solving  nonlinear 
equation 

SL  =  search  length  for  locating  first  change  of  sign 
in  nonlinear  equation  solver  (subroutine 
ROOMID) 

CONST  =  constant  in  Eqn  (1) 


ZOMR  =  ^ 


Z  = 


USTAR  =  1  -  R* 

In  general,  the  main  program  reads  in  the  necessary 
data,  computes  using  ROOMID  and  ZORC, 


then  uses  ROOMID  and  RTF  to  calculate 


11 
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FUNCTION  RTF(2) 

Significant  Variable  names: 


SCRIT  =  ^ 


H  =  interval  size 


RTF  is  the  function  described  by  Eqn.  (1)  (a  function 


).  The  tail  of  the  integral 


of 


computed  first,  using  QUAD  until  the  contribution  to  the 
integral  becomes  insignificant.  GLAG  is  then  used  to 
squeeze  the  last  drop,  if  any,  out  of  the  tail.  The  inte- 
gral  for  is  then  computed  using  QUAD. 

C.  FUNCTION  F  (S) 

F  computes  the  integrand  in  Eqn.  (1). 

D.  FUNCTION  ZORC  (X) 


ZORC  (X)  =  ^ (x) 


+  Rq  -  1 . 


ZORC  is  used  to  calculate 


E.  QUAD  (A,  B,  AR) 


QUAD  evaluates  the  integral  of  F  from  A  to  B 
using  three  point  Gauss-Legendre  quadrature. 


F.  GLAG  (A,  AR) 

GLAG  evaluates  the  integral  of  F  from  A  to  infi¬ 
nity  using  five  point  Gauss-Laguerre  quadrature. 

G.  CON  (C,  N,  CONST) 

CON  evaluates  the  constant  term  in  Eqn.  (1). 


12 


H.  FUNCTION  DN  (X) 

DN  is  the  normal  distribution  function  ^  , 

I.  FUNCTION  PINT  (X) 

PINT  calculates  for  large  negative  arguments, 

'J*  ROOHID  (G,  X.  SL,  FRE ) 

ROOHID  solves  the  equation:  G  (X)  =  0  using 
the  midpoint  method.  X  is  the  initial  guess  at  the  root; 

SL  is  the  search  length.  ROOHID  assumes  that  SL  will 
point  in  the  true  direction  of  the  root;  i.e.,  if  X  initi¬ 
ally  overestimates  the  true  root,  SL  must  be  negative. 
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APPENDIX 


Use  of  Tables: 

The  entries  in  the  body  of  the  tables  give 

A 

for  a  given  sample  size  n  and  confidenced  reliability  R^. 
Eqns.  (3)  and  (5)  are  used  to  calculate  normal  and  log¬ 
normal  mission  life  respectively.  If  mission  life  is  the 
given  quantity  and  confidenced  reliability  is  desired,  one 
could  use  Lagrangian  interpolation,  but  the  simplest  thing 
to  do  would  be  to  draw  a  graph  of  R^.  versus  ^ 
desired  n.  An  example  is  given  for  n  =  2* 

A  graph  of  versus  n  is  included  for  R^  =  .999 

and  C  =  .90.  This  graph  indicates  that  a  considerably 
rapid  increase  in  information  is  obtained  by  increasing  n 
from  two  to  five,  while  increasing  n  beyond  10  produces 
information  at  a  very  slow  rate. 

Finally,  it  must  be  noted  that  maximim  likelihood 
(not  unbiased)  parameter  estimates  are  to  be  used  with  these 
tab les . 
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C  =  90% 


2 

3 

4 

5 

6 

7 

8 

9 

10 

15 

20 

25 

999 

34.764 

11. 820 

8.232 

6.  833 

6.086 

5.619 

5.297 

5 . 060 

4.879 

4.  363 

4.  113 

3.962 

995 

28.972 

9.911 

6.914 

5.742 

5.115 

4.722 

4.451 

4.252 

4.099 

3.663 

3.451 

3.322 

99 

26. 163 

8.990 

6.280 

5.217 

4.647 

4.290 

4.044 

3.862 

3.723 

3. 325 

3.131 

3.013 

975 

22.042 

7.647 

5.354 

4.451 

3.966 

3.661 

3.450 

3.294 

3.  174 

2.831 

2.664 

2.562 

95 

18.512 

6 . 505 

4.569 

3.801 

3.387 

3. 126 

2.944 

2.811 

2.707 

2. 411 

2.265 

2. 176 

925 

16.232 

5.772 

4.064 

3.383 

3.015 

2.781 

2.619 

2.499 

2.407 

2.139 

2.008 

1.927 

90 

14.500 

5.215 

3.681 

3.066 

2.732 

2.520 

2.372 

2.262 

2.177 

1.932 

1.  811 

1.737 

875 

13.083 

4.760 

3.367 

2.806 

2.500 

2.305 

2.16  8 

2.067 

1.989 

1.  762 

1.649 

1.579 

85 

11. 877 

4.371 

3.099 

2.583 

2.300 

2. 120 

1.994 

1.900 

1.  827 

1.615 

1.509 

1.  444 

30 

3. 859 

3.235 

2.933 

2.492 
2. 115 

1.  872 
1.685 
1.532 
1.399 
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2 

3 

4 

5 

6 

7 

999 

69.687 

16.971 

10.64 

8.387 

7.243 

6.548 

6. 

995 

58.084 

14.241 

8.947 

7.058 

6.097 

5. 512 

5. 

99 

52.458 

12.924 

8.132 

6. 419 

5 .545 

5.014 

4. 

975 

44.205 

11.006 

6.945 

5.488 

4.  742 

4.  288 

3. 

95 

37. 137 

9.377 

5.9  40 

4.699 

4.062 

3.672 

3. 

925 

32 .573 

8.  331 

5.295 

4.  193 

3.625 

3.276 

3. 

90 

29.107 

7.539 

4.  806 

3.809 

3.293 

2.9  76 

2. 

875 

26.273 

6.891 

4.406 

3.495 

3 .022 

2. 730 

2. 

85 

23. 860 

6.338 

4.064 

3.226 

2 .789 

2. 519 

2. 

C=95% 
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10 

15 

20 

25 

30 

5.742 

5 . 485 

4.  769 

4.430 

4.228 

4.091 

4.832 

4.615 

4.010 

3.722 

3.550 

3. 434 

4.394 

4.196 

3.644 

3.381 

3.223 

3.116 

3.757 

3.586 

3.110 

2. 883 

2.7  46 

2.653 

3.215 

3.068 

2.656 

2.458 

2.339 

2.258 

2.867 

2.735 

2.364 

2.185 

2.076 

2.00  3 

2.603 

2.482 

2.141 

1.976 

1.  876 

1.  808 

2.386 

2.274 

1.958 

1.  804 

1.  711 

1.647 

2.199 

2.096 

1.  800 

1.657 

1.  569 

1.509 

8 

080 

118 

654 

980 

407 

040 

760 

531 

335 


t 


I 


\ 

1 


i 


17 


35 


30 


/fc  = 

C  =  -lo 


25 


5 


*  *  ■»  4  »  * • 1— ♦ * t  < — 

10  /5  20  25  30 

/n 
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100  DIMENSION  R(50) 

110  COMMON  RC>C0NST.N»RN^ZJMR»SQN»RN02 

120  COMMON  NM2*ZZ#RINT 

130  COMMON  S2PI#NINT( 10) > IT 

140  EXTERNAL  ZORC»RTF 

150  CALL  OPENF(  1  ,  ••TBL4'*) 

160  READC 1 , ) ,N,NR,NIT 
170  READC 1 » ) » CR( I ) , 1=1 ,NR) 

180  READC 1 » ) . (NINTC I ) » I =1 >NI T ) 

190  READC 1 , ) , C. FRE,SL 
200  PRINT  500. N,C 

210  500  FORMATC  1  4HN0  .  OF  PO  INTS  =  .  1  3// 1  7riC0NFI  DENCc.  LE\/c.L=  /  ^  5 . 3/ / ) 
220  CALL  CJNCC.N. CONST) 

230  PRINT  600. CONST 

240  600  FORMATC 6HC0NST=.E16.8/) 

250  RN=N 

260  SON=SQRTCRN) 

270  PI=3. 1415927 
280  S2PI=S0RTC2.*PI ) 

290  RND2=RN/2. 

300  NM2=N-2 
310  DO  1  IR=1.NR 
320  SLL=SL 
330  RC=RCIR) 

340  Z0tyR=0. 

350  CALL  R00MIDCZ0RC.Z0MR.-2..FRE) 

360  PRINT  100. RC 

370  100  FORMATC24HCONFIDENCED  RELI ABI LI TY=. F7 . 5/ ) 

380  PRINT  200.ZOMR 

390  200  F0RMATC8HZC 1-RC)=.E16.8/) 

400  Z=ZOMR 

410  DO  1  I1T=1.NIT 

^20  RINT  =  N1NTC I  IT) 

430  IT=1IT 

4/<0  CALL  ROOMIDCRTF.Z.SLL.FRE) 

450  ZSTAH=-Z 

460  PRINT  300.ZSTAH 

470  300  F0RMATC9HZCRSTAR)=.E16.8/) 

480  USTAR=DNCZ) 

490  PRINT  400.USTAR 

500  400  FORMATC 6HUSTAR=. El  6.8/) 

510  SLL=.01*Z 
520  Z=Z-SLL 
530  1  CONTINUE 
540  END 
550C 
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570  FUN'CTIOiNJ  RTF(Z) 

530  COP'/, ..IN  RC>CJNST»N,ki\%Z)Mh>SGN>f<NJ^i 

590  CO/iMON  iM/2>ZZ^  RINT 

600  COP'MON  S2PI  >NINT(  10)  >  I  T 

610  ZZ=Z 

620  SCRIT=ZOMR/Z 
630  AREA1=0. 

640  H=SCRIT/CRINT-.5) 

650  A=SCRIT+.5*H 

660  1  B=A+H 

670  CALL  QUAD(A»B^AR) 

680  AREA2=AREA1+AR 

690  IF(AREA2-AREA1 )3>3#2 

700  2  AREA1=AREA2 

710  A=B 

720  (-0  TO  1 

730  3  CALL  (LA(  <A>AR) 

7/10  AREA  =  AREA1 +AK 
750  NNT=NINTCIT) 

760  DO  4  I=1>NNT 
770  RI=1 
780  A=CRI -  1  . ) *H 
790  B=A-'-H 

800  CALL  QOAD<A,B/AR) 

810  AREA=AREA+AR 

820  ^  CONTINUE 

830  RTF=AREA-CJNST 

840  RETURN 

850  END 

B60C 

870C 

880  FUNCTION  FCS) 

890  COMMON  fiC» CONST* N > RN, ZOMR, SON » RN 02 
900  COMMON  NM2*ZZ>RINT 
910  ARGPH= SON* ( S*ZZ-Z OMR) 

920  PHI=DN(  AR(.PH) 

930  F=PH1 *S**NM2*EXP( -RN02*S*S) 

940  RETURN 
950  END 
960C 
970C 

980  FUNCTION  ZORCCX) 

99  0  COP'MXN  RC 

1000  ZORC=DN(X)+RC- 1 . 

1010  RETURN 
1020  END 
1030C 
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1040C 

1050  SUBROUTINE  OUAD<A,B,AR) 

1060  BMA2=(B-A)/2. 

1070  Xl=. 77459667 
1080  A1=BMA2*<X1  +  1 .  )  +  A 
1090  A2=BMA2*C -XI + 1 . ) +A 
1 100  A0=(B+A)/2. 

1 1 10  G1=F(A1 ) 

1120  G2=FCA2) 

1 1 30  GO=FCAO) 

I  140  AR=BMA2*(5.*Gl+8.*G0+5**G2)/9. 

1150  RETURN 

1160  END 

I I  700 
1  1800 

1190  SUBROUTINE  GLAGCA,AR> 

1200  DIMENSION  X(5)^WEX(5) 

12100  PROGRAM  OOMPUTES  INTEGRAL  FROM  A  TO  INFINITY  OF  F(X) 

12200  USING  5  POINT  GAU SS-LAGUERRE  QUADRATURE 

1230  X( 1 )= .26356032 

1240  X(2)=l .4134031 

1250  X(3)=3. 5964258 

1260  X(4)=7 .0858100 

1270  X(5)=12. 640801 

1280  WEX( 1 )= .67909404 

1290  WEXC2)=1 . 6384879 

1300  WEXC3)=2. 7694432 

1310  V.'EX(4)=4. 31  56569 

1320  WEX( 5)=7 .2191864 

1330  AR=0. 

1340  DO  1  1=1,5 

1350  A1=X( I )+A 
1360  G=F(A1 ) 

1370  AR  =  AR-^G*WEX(  I  ) 

1380  1  CONTINUE 
1390  RETURN 
1400  END 
1410C 
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1430  SUBROUTINE  CON ( N . CONST ) 

1440C  PROGRAM  COMPUTES  < 1 -C ) GAK< C N- 1 ) /2 ) / ( 2< N/2 ) * * ( ( N- I ) /2 ) ) 
1450  RN=N 
1460  N2=N/2 
1470  RN02=RN/2. 

1480  RN2=N2 

1490  IF<RN02-RN2> 3^ 1 ^ 3 
1500  1  PI=3. 1415927 
1510  P=1 . 

1520  D'.l  2  I  =  WN2 
1530  RI=I 

1540  P=P*(2.*RI-1.) 

1550  2  CONTINUE 

1560  CONST=P/( (2. ♦RN2- 1 . ) * ( 2 . ♦RN2 ) * ♦N2 ) * ( 1 . -C ) *  SORT ( P I  * KN2 ) 
1570  RETURN 
1580  3  P= 1 . 

1590  DO  4  1=1 /N2 
1600  RI=I 
1610  P=P*RI 
1620  4  CONTINUE 

16  30  C0NST=P*2.**(N2-  1  )  /  f  RN2*  C 2 .  *RN2+  1  .  )**N2)*  (  1.-0 

1640  RETURN 

1650  END 

1660C 

1670C 

1680  FUNCTION  DN<X) 

1690  IFCX+5. ) W 1>2 
1700  1  DN=PINT(X) 

1710  RETURN 

1720  2  IF(X-6. ) 4> 3^ 3 

1730  3  DN=1 . 

1740  RETURN 
1750  ^  AX=ABS(X) 

1760  S2PI=2. 5066283 
1770  P=. 2316419 
1780  Bl=. 31938153 
1790  B2=- .35656378 
1800  B3=l. 7814779 
1810  B4=-l .8212560 
1820  B5=l .3302744 
1830  T= 1 ./ C 1 . +P*AX) 

1840  S=T*CB1+T*(B2+T*CB3+T*(B4+T*B5) ) ) ) 

1850  DN=S/S2PI*EXP(-.5*X*X) 

1860  IFCX)5»5>6 
1870  5  RETURN 
1880  6  DN= 1 . -DN 
1890  RETURN 
1900  END 
19  IOC 
1920C 
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1930  FUNCTION  PINT(A) 

1940  A2=A*A 
1950  A3=A*A2 
1960  Sl=< 1 . -A2)/A3 
1970  S2=S1 
1980  AI  1=  1  ./A3 
1990  DO  1  I=2>200 
2000  RI=I 

2010  AI2=-(2.*RI-*- 1  .  )*AI  1/A2 

2020  IFCABSCAI2)“ABSCAI 1))2j4,4 

2030  2  S2=S2+AI2 

2040  IFCSl -S2)3. 4,3 

2050  3  S1=S2 

2060  Al 1=AI2 

2070  1  CONTINUE 

2080  4  C=2. 5066283 

2090  PINT=EXP(-A2/2. >♦32/0 

2100  RETURN 

21  10  END 

2120C 

21  30C 

2140  SUBROUTINE  ROOMI D( G, X, SL , FRE ) 
2150  EXTERNAL  G 
2160  X1=X 
2170  F1=G(X1) 

2180  3  X2=X1+SL 
2190  F2=GCX2) 

2200  IF(F1*F2)5,5,2 

2210  2  X1=X2 

2220  F1=F2 

2230  GO  TO  3 

2240  5  X=.5*(X1+X2) 

2250  FX=6CX) 

2260  IF(FX*F1 )7, 1 3, 6 

2270  6  XI =X 

2280  F1=FX 

2290  GO  TO  10 

2300  7  X2=X 

2310  10  A1=ABSCX1 ) 

2320  A2=ABSCX2) 

2330  IF(A1-A2)8,B,9 
2340  8  AMIN=A1 
2350  GO  TO  1 1 
2360  9  AMIN=A2 

2370  11  PE=.5*ABS(X2-X1)/AMIN 
2380  IF(PE-FRE) 12, 12, 5 
2390  12  X=.5*(X1+X2) 

2400  13  FX=F(X) 

2410  RETURN 
2420  END 


23 
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